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ABSTRACT 

The discovery of radio pulsars in compact orbits around Sgr A* would allow an unprecedented 
and detailed investigation of the spacetime of the super massive black hole. This paper shows that 
pulsar timing, including that of a single pulsar, has the potential to provide novel tests of general 
relativity, in particular its cosmic censorship conjecture and no-hair theorem for rotating black holes. 
These experiments can be performed by timing observations with 100 /is precision, achievable with 
the Square Kilometre Array for a normal pulsar at frequency above 15 GHz. Based on the standard 
pulsar timing technique, we develop a method that allows the determination of the mass, spin, and 
quadrupole moment of Sgr A*, and provides a consistent covariance analysis of the measurement 
errors. Furthermore, we test this method in detailed mock data simulations. It seems likely that only 
for orbital periods below ~ 0.3 yr is there the possibility of having negligible external perturbations. 
For such orbits we expect a ~ 10~^ test of the frame dragging and a ~ 10~^ test of the no-hair theorem 
within 5 years, if Sgr A* is spinning rapidly. Our method is also capable of identifying perturbations 
caused by distributed mass around Sgr A*, thus providing high confidence in these gravity tests. Our 
analysis is not affected by uncertainties in our knowledge of the distance to the Galactic center, Rq. 
A combination of pulsar timing with the astrometric results of stellar orbits would greatly improve 
the measurement precision of Rq. 

Subject headings: black hole physics; Galaxy: center; pulsars: general 



1. INTRODUCTION 

^ ■ One of the most intriguing results of general relativ- 
ity (GR) is the uniqueness theorem for the stationary 
black hole s olutions of the Einstein-Maxwell equations 
^-H (see Heuslerl [19981 . and references therein). This unique- 
ness theorem states that (under certain conditions) all 
stationary electrovacQ black hole spacetimes with a non- 
degenerate horizon are described by the Kerr-Newman 
metric. It implies that in GR all stationary black holes 
are parameterized by only three parameters: mass (M), 
spin (S) and electric charge ("black holes have no hair"). 
. ^ _ All uncharged black hole solutions are described by the 
k> . Kerr metric and, therefore, uniquely determined by M 
^ and S. Astrophysical black holes are believed to be the 
d result of a gravitational collapse. During this collapse 
all the properties of the progenitor, apart from mass and 
spin, are radiated away by gravitational radiation while 
the gravitational field asymptotical ly approaches its sta- 
tionary configuration ()Pricel[T972al lbl). The outer space- 
time of an astrophysical black hole should therefore be 
described by the Kerr metric^ Since the Kerr metric 
has a maximum spin at which it still exhibits an event 
horizon, Pen rose's cosmi c censorship conjecture (CCC) 
within GR (;Penrosell979f ) requires the dimensionless spin 



^ Electrovac spacetimes are the solutions of the Einstein-Maxwell 
equations. 

^ Strictly speaking, this is only true for a certain approximation 
since, to some extent, astrophysical black holes will be influenced 
by nearby masses (accretion, orbiting objects). We will address this 
issue for Sgr A* in this paper at the end of the discussion Section. 



parameter x to satisfy 
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A measured value for x that exceeds 1 would pose a 
serious problem for our understanding of spacetime, since 
this would indicate that either GR is wrong or that a 
region may be visible to the outside universe, where our 
present understanding of gravity and spacetime breaks 
down. 

As a result of the no-hair theorem, all higher multipole 
moments (/ > 2) of the gravitational field of an astro- 
physical blac^hole_can be expressed as a function of M 
and S (H anseni [1971 ) . In particular, the dimensionless 
quadrupole moment q satisfies the relation 
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A measurement of the quadrupole moment, in combina- 
tion with a mass and a spin measurement, would there- 
fore provide a test of the no-hair theorem for Kerr black 
holes. 

Some of the clearest evidence for the existence of black 
holes comes from the mo nitoring of ^ 30 stellar orbits in 
the center of our Gala xv (iSchodel et al.ll200l iGhez eFall 
120081 : iGillessen etHI |2009D, where the shortest orbital 
period, Pb, is 16 years. Known by its radio nomencla- 
ture of Sgr A*, current estimates put the mass of this 
black hole to around 4 x 10^ M©. A black hole of that 
size at a distance of 8 kpc is an ideal laboratory for black- 
hole physics, strong-field gravity and, in particular, a test 
of the no-hair theorem for Kerr black holes (e.g. .Psaltisi 
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2QQ8I : iJohannsen fc Psaltisll2Qll[ ). It has been shown by 
Willi ()2QQ8f ) that the discovery of stars in highly eccentric 
(e ~ 0.9) orbits very close to Sgr A* (Pb < 0.1 yr) could 
be used to test the general relativistic no-hair theorem. 
This experiment requires an astrometric precision at the 
level of 10 /ias, which seems achievable with the upcom- 
ing infrared- astrometry experiments, such as GRAVITY 
(jEisenhauer et al1l2009l ). At a distance of 8 kpc, an angle 
of 10 //as corresponds to a length scale of ~ 10 km. 

On the other hand, if close-in pulsars could be found 
and tracked in their orbits, even for those with poor tim- 
ing precision the time-of- arrival (TOA) for a (integrated) 
pulsar signal can be measured with an uncertainty of 
a few milliseconds, corresponding to a light-travel dis- 
tance of ~ 10^ km. Moreover, a phase-connected solu- 
tion with an appropriate timing model allows a deter- 
mination of the pulsar orbit, w hich provides eve n more 
pre cision. Hence, as sho wn by IWex fc KopeikinI () 19991 ) 
and lKramer et al.l (|2004f ). a pulsar in orbit around the su- 
permassive black hole in the Galactic center (GC) would 
provide an ideal probe to measure the mass, the spin, and 
the quadrupole moment of Sgr A*, and consequently test 
the no-hair theorem for Kerr black holes. In their discus - 
sion on "bumpy black holes" , IVigeland fc HughesI (j2010i) 
suggest that a pulsar in orbit around a black hole could 
be used for mapping its multipole moment structure. In 
a recent publication 'Angelil et al.' (^2010) have discussed 
the importance of frame dragging and the quadrupole 
moment for stars and pulsars in orbit around Sgr A*, 
based on numerical integration of geodesies in a Kerr 
spacetime. For pulsar s, however, the resul ts in this paper 
are only indicative Q I Wang etaH ()2009al lbl) have shown 
that if a pulsar is found in a very close orbit around 
Sgr A* (e.g. Pb ^ Iday), observers at the Earth can 
receive additional pulses travelling along a path that is 
strongly bent by the gravitational potential of the black 
hole. The exploitation of this information would provide 
unique constraints on the strongly curved spacetime ge- 
ometry near Sgr A*. 

In this paper we will demonstrate how the mass, the 
spin (magnitude and orientation) and the quadrupole 
moment of Sgr A* can be determined from timing a pul- 
sar in a sufficiently tight orbit around the supermassive 
black hole. Our analysis is based on simulated TOAs 
and a timing model that allows for a phase-connected 
solution, consistently accounting for the relativistic ef- 
fects in the motion of the pulsar and the propagation of 
the radio signals. Based on this, we can determine the 
expected precision for the individual parameters while 
accounting for all the correlations between the parame- 
ters. In all discussions and simulations we focus on the 
leading order in the individual effects of interest. We are 
well aware that in an actual timing model for orbits close 
to Sgr A* we need to account for higher order effects in 



^ The "Angelil et al." ("2010") results for a pulsar are not derived 
from a consistent covariance analysis based on a timing model that 
incorporates all the relevant effects simultaneously. Moreover, pul- 
sar timing is treated as a radial velocity measurement experiment, 
which is incorrect. In fact, pulsar timing makes use of phase- 
connected solutions for the rotational phase of a pulsar leading to 
a precision in the parameter estimations that can be several orders 
of magnitude better (Lorimer & Kramer 2005). Also, the timing 
precision assumed by Angelil et al. ( 2010) seems too optimistic for 
a GC pulsar, as will become clear from the analysis presented in 
this paper. 



the orbital motion and signal propagation. For many of 
the effects discussed belo w, higher order terrns have al- 
ready been calculated fe.g. Damour fc S chafer" 1 988': ' We:?3 
1995; Kopeikin 1997; Konigsdorffer fc Gopakumar 20^ 
However, to estimate the expected precision and covari- 
ances in the parameter determination, it is sufficient to 
use a timing model that combines just the leading terms 
of all contributions relevant here. The second order terms 
contribute at a ~ /^q level, w here 0o is the orbita l veloc - 
ity parameter introduced by iDamour fc Taylorl ()1992l ). 
For a test particle in orbit around a mass M 
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As an example, for an orbital period of 0.1 yr we find 
pQ ^ 10~^. At this point it is worth mentioning, that 
for the most relativistic bin ary pulsar, the double pulsar, 
one finds 4.3 x 10"^ (|Kramer fc Wexl[2009h . which 

nicely illustrates how much more relativistic a pulsar in 
a Pb ^ 1 yr orbit around Sgr A* would be. 

Previous studies suggest that about 1000 pulsars can 
be expected t o be orbiting Sgr A* with periods less 
than 100 years (ICordes fc Laziolll997H(^nzel et, al.ll2nn3l : 
IPfahl fc Loebl[2nnl FI. and some of them may be associ- 
ated with remnants of the observed S-star population 
in the neighborhood. Deep pulsar searches towards the 
GC region were already conducted with a few radio tele- 
scopes (Effelsberg, Green Bank, Parke s ) at frequencies 
up to 15 GHz (e.g. iKramer et all 120001 : iMacquart et al.l 
"2010). Five pulsars were found no more than 200 pc away 
from Sgr A* (Johnston et al. 2006; Deneva et al. 20i)9, 
Kramer et al. private communication), which is consis- 
tent with the estimated large pulsar population within 
that region. However, they are not close enough to the 
supermassive black hole to probe its gravitational field. 
In this paper we will focus on timing observations of such 
pulsars, and show how far they could take us in probing 
the gravitational field of Sgr A*, provided the system is 
found to be sufficiently free of external perturbations. 
The main purpose therefore is the development of the 
methodology, and the estimation of its potentials in test- 
ing the Kerr nature of Sgr A* based on mock data sim- 
ulations. For further elaboration on either the existence 
of pulsars in close orbits around Sgr A* or the search for 
them , we refer to the rich literature (jLorimer fc Kramerl 
l2005l ) and future work in progress. In Section [2] we dis- 
cuss the expected timing precision and the orbital peri- 
ods required for our measurements. Section [3] presents 
the various relativistic effects that can be used to de- 
termine the mass of Sgr A* and, based on simulations, 
the expected precision in the mass measurement. In Sec- 
tion [4] we show how the spin can be extracted from the 
timing measurement, how this information can be used 
to test the CCC, and how the presence of a distributed 
mass in the vicinity of Sgr A* would affect this measure- 
ment. Section [5] provides the details on the quadrupole 
measurement and the no-hair theorem test. In Section [6] 



The result from the GC survey bv IMacquart et al.l (|20Tol ) in- 
dicated that the actual number of such pulsars may be less. 
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we summarize and discuss our findings. 

2. TIMING A PULSAR IN ORBIT AROUND Sgr A* 

Pulsar timing involves measuring the TOAs of a pul- 
sar's pulses aiidm them on a timescale of years 
fe.g. lTavloi][T992h . The precision of TOA measurements 
of young pulsars near the GC, by future telescopes, will 
mainly be limited by three effects that have significantly 
different dependencies on observing frequency: firstly the 
signal-to-noise ratio of the measured pulses, secondly the 
pulse phase jitter intrinsic to the pulsar, and thirdly the 
changes in pulse shape caused by interstell ar scintillation 
(jCordes Shannonl l20inl : iLiu et alJ[20TTh . The first ef- 
fect is independent of frequency under our assumption 
that the pulse width does not change with frequency. In 
reality, the pulse width does evolve but that is secondary 
to the overall timing error. The second and third effects 
are strongly frequency dependent, due to the steep pulsar 
spectrum and the pulse broadening caused by the large 
amount of scattering from the high electron density in 
the ionized gas near the GC. The strong dependence of 
scattering on frequency (oc /~^, see the next paragraph) 
implies that observations need to be made at much higher 
frequencies than are typically used for pulsar timing. 

There have been previous studies on optimizing the ob- 
servational frequency for the purp ose of pulsar searches 
towards the central parsec region (jCordes fc Laz io'^TOQ?: 
iMacquart et al. 2010). Fig.[T]shows the estimated timing 
precision for a canonical pulsar near Sgr A* as a function 
of the observing frequency. The calculation of the achiev- 
able TOA uncertainty (Jtoa accounts for three contribu- 
tions: 

2 2,2,2 / A\ 

^TOA = ^rn + + ^scint ' iV 

Here am, crj, and cTgcint represent the uncertainties con- 
tributed by radiometer noise, pulse phase jitter and inter- 
stellar scintillati on, respectively, whi c h can be calculated 
by following e.g. iCordes fc Shannon! ()2010f ). Specifically, 
we use a spin period P = 0.5 s, an intrinsic pulse width 
Wi = 10 ms, and a period-averaged flux density 51400 = 
1 mJy at 1.4 GHz. For a 100-m diameter dish and the 
SKA, we use a gain of 2K Jy~^, 100 K Jy~^, and an in- 
tegration time of 4hr, 1 hr, respectively. Two different 
spectral indices of the pulsar flux de nsity, which t ypifies 
many of these measured for pulsars (jMaron et al.| [2000). 
are used in our calculations. The scattering time scale is 
estimated to be Tgcat ^ 2.3 x 10^ ms at 1 GHz, as derived 
from the observed scattering diameter of Sgr A* and the 
estimated location of the scattering material along the 
line of sight, the latter as incor porated in the NE20Q 1 
model = 6 = and D=8.5 kpc, ICordes fc La^l2002f ). 
For this large a mount of scattering, we use a scaling of 
''"scat cx: (e.g. Loh mer et al.|[2QQlh rather than the of- 
ten used Kolmogorov scaling Tgcat cx: /"^"^ (e.g. iRickettI 
[1990) , because the dominant length scale is less than the 
inner scale of the wavenumber spectrum for the electron 
density. Note that all potential pulsars with close orbits 
of interest for the GR tests will be seen along essentially 
the same line of sight as Sgr A*, so one can assume that 
their lines of sight will have the same scattering char- 
acteristics. The system temperature (e.g. :^ 40 K at 
15 GHz) is calculated by summing the radio background, 
receiver temperature and emission of the atmosphere. It 
clearly follows from Fig.[Tl that with a radio telescope like 



SKA TOA uncertainties of below 100 /is seem likely for 
an observational frequency above 15 GHz, similar to the 
result of optimized searching frequency. A detection of 
MSPs in the Galactic centre is unlikely (I Cordes fc Laziol 
119971 : LMacquart et al.' '2010) , so they are not considered 
in the following simulations. However, we will show that 
the black hole properties can already be extracted by 
finding and timing a relatively slow pulsar. If a MSP 
were to be found after all, the experiment may reach a 
correspondingly higher precision. 
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Fig. 1. — Predicted TOA measurement precision of a young pul- 
sar near the GC for two different spectral indices a. The curves 
take into account pulse phase jitter intrinsic to the pulsar, pulse 
broadening from scattering along the entire line of sight, and from 
the finite number of scintles included in the measurement. Scatter- 
ing is dominated by a region of high plasma density that surrounds 
the GC. We assumed a four hour integration time using a 100-meter 
radio telescope and a one hour integration time using an SKA-like 
telescope, both with a highest operating frequency of 30 GHz and 
a bandwidth of 1 GHz. It is found that observational frequen- 
cies above 15 GHz favors pulsar timing observation, where 100 /iS 
TOA precision seems achievable, in particular with the SKA. The 
parameters used for this calculations can be found in the text. 

Precision of long-term timing of a young pulsar is often 
limited by irregularities of the pulsar's rotation, suppos- 
edly associated with either the in ternal super-fluid flux 
(e.g. 'Mela tos fc Warszawski l'2QQ9'). or external magnetic 
field activity (e.g. iLvne et " al. 2010). The amplitude of 
the low-frequency noise resulting from these instabilities 
varies from about 10 /is to 100 ms (Hobbs et al. 2010), 
and in some cases the noise can be model ed by follow- 
ing the approach proposed by .Lvne et all ()2010f ) to im- 
prove the timing precision by order of magnitudes. Con- 
sequently, a TOA precision of 100 jas is a reasonable fidu- 
cial value, which we will use in our simulations below. 

Although the purpose of this paper is to discuss poten- 
tial gravity tests with a pulsar in orbit around Sgr A*, 
provided the system is found to be sufficiently clean, we 
nevertheless would like to complete this Section with 
a brief discussion on possible effects t hat cou ld com- 
plicate or eve n spoil t hese tests. Merritt et al.j (j2010l) 
and Sadeghia n fc Will (2011) have shown that for or- 
bits with an orbital period Pb larger than 0.1 yr, it be- 
comes likely that the distribution of stars in the vicinity 
causes "external" perturbations of the orbital motion of 
the pulsar and prevent a clean test of the no-hair the- 
orem or even a measurement of the Lense-Thirring ef- 
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feet. In order to evaluate the si gnificance of the per - 
turbation, following the analysis of lMerritt et al.l ()2Q1Q| ). 
in Fig. [2] the relation of precessional timescale against 
orbital size is presented for four different contributions: 
the pericenter advance, the frame dragging effect, the 
black- hole quadrupole, and the surrounding mass distri- 
bu tion. Here we assum e 10^ (the highest number applied 
m iMerritt et aIll2QlQt ) one solar mass objects isotropi- 
cally distributed wit hin 1 mpc around Sgr A*. Similar to 
IMerritt et al.l ()2Q1Q| ) we do not consider the influence of 
objects outside the central 1 mpc region. One can see 
that for wide orbits (Pb ~ 10 yr) the pericenter advance 
is still significantly larger than the precession caused by 
external perturbations. This suggests that for orbital 
periods less than about 10 years the measured uj can be 
used to tightly constrain the black hole mass. The frame 
dragging will be dominant over the stellar noise if the 
orbital period is less than 0.5 years, while only for or- 
bital periods < 0.1 yr the (secular) contribution of the 
quadrupole moment is expected to be significantly above 
the external perturbation. We note that the assump- 
tions applied to calculate the precessional timescale by 
stellar perturbation may not be secure, as the actual stel- 
lar components and mass distributions within the cen- 
tral pc (especially the central mpc) are still not fully 
understood. For instance, the precessional torque could 
be larger if there exists a high fraction of massiv e ob- 
jects near Sgr A* due to mass segregation CO Tearv et al.l 
|2009; Keshet et al.ll2009l : IKocsis fc Tremain e 2011), or a 
significant anisotropy in the distribution of the surround- 
ing masses. In fact, it has been argued bv IMerritt et al.l 
(j2010.) that a high fraction of 10 Mq black holes in this 
region would make astrometric tests of general relativity 
problematic at all radii. Concerning pulsars, however, 
as will be discussed in Sections 14.31 and 15.21 the grav- 
itomagnetic and quadrupolar field of Sgr A* will result 
in unique features in the timing residuals, which can be 
tracked well with high precision timing observations, and 
one can still expect to be able to extract the Sgr A* spin 
and quadrupole moment from the timing data, to some 
extent. This however depends strongly on the details 
of the external perturbations, which only will be known 
once a pulsar is discovered in that region. 

3. MASS MEASUREMENT 

The current best estimate for the mass of Sgr A* gives 
4.30 ± 0.20(stat) ± 0.30(sys) X 10^ (Gillessen et al. 
l2009f ]Fl. The proposed method has the potential to im- 
prove the measurement accuracy by a factor of ^ 10^. 
This is possible pulsar timing. Indeed, the most precise 
mass measurements for stars (other than the Sun) come 
from pulsar tim ing observations (|Jacobv et al.l |2005l; 
iFreire et al.ll2QTTI ). Those are achieved in binary pulsar 
systems, where in addition to the Kepler ian parameters 
one can determine a set of post-Keplerian (PK) param- 
eters as theory independent relativistic corrections. In 
any theory of gravity the PK parameters are functions 
of the two a priori unknown masses of the system, which 
can be deter mined once two PK para meters have been 
obtained (see iDamour" fc Tavlor|[T992l for definitions of 
the Keplerian and PK parameters). Since in our case 

^ The main uncertainty is from the hmited knowledge of the 
distance to the GC. 
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Fig. 2. — Timescales of secular orbital precession for a pulsar in 
orbit around Sgr A* as a function of orbital period Pb (semi-major 
axis a). The letters M, S, Q and P stand for the contribution 
by the mass monopole (pericenter advance), the spin (frame drag- 
ging), the quadrupole moment, and stellar perturbation, respec- 
tively. Here we assume an orbital eccentricity of 0.5 and 10^ one 
solar mass objects within 1 mpc around Sgr A*. As a comparison, 
the Schwarzschild radius of Sgr A* is ~ 4 x 10"'^ mpc. 



the mass of the pulsar can be neglected in comparison 
to the mass of the black hole, in general one PK param- 
eter is sufficient to estimate the mass of Sgr A* with a 
precision at the 10~^ level. In the following we briefly 
discuss three relativistic effects that can be used for a 
mass determination, and present the results of extensive 
simulations at the end of this Section. As mentioned in 
Section [H in a discussion of measurement precision it is 
sufficient to keep the first order terms in the description 
of these effects. 

3.1. Post-Keplerian parameters and mass determination 

In eccentric binary pulsars the precession of periastron^ 
cj, is usually the first PK parameter that can be measured 
with high precision. It allows the determination of the 
total mass of the system, which in our case can be equaled 
with the mass of the black hole Mbh- To first order one 
finds (^Robertson ■ 1 938i ) : 
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(5) 

where e is the orbital eccentricity. As an example, for 
a O.Syr orbit with an eccentricity of 0.5 the orbit pre- 
cesses with a rate of about 2.7 degrees per year. After 
five years of weekly observations with a timing uncer- 
tainty of 100 /is, this precession will be measured with a 
fractional precision of better than 10~^. This, however, 
is not the precision with which the mass of Sgr A* can 
be determined. If the black hole is rotating, a significant 
fraction of the p ericenter precession can co me from the 
frame dragging ([Barker fc 0'Connelllll975[ ). Depending 
on the orientation and the spin of the black hole, this 
could be up to about 1% of the total precession. A mea- 
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surement of the spin and the orientation of the black 
hole would allow to correct for this Lense-Thirring con- 
tribution co'LT- But as we will show later, we will use 
the observed oj and the mass measurement from other 
relativistic effects to calculate cOlt and use it in the spin 
determination {uj = uju + ^lt)- 

The Einstein delay is a combination of the second order 
Doppler effect and gravitational redshift. From its am- 
plitude 7e, which is a PK parameter, one can determine 
the mass of the black hole, since (to first post-Newtonian 
order) (£landford & Teukolsky 1976) 



7e ^2e 
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For a 0.3 yr orbit with an eccentricity of 0.5 the ampli- 
tude of the Einstein delay will be of order 800 seconds. 
However, the Einstein delay is a priori not separable from 
the Roemer dela}0, and can only be measured with suffi- 
cient accuracy after some time, when the relativistic pre- 
cession has changed the orbital orientation sufficiently. 
For a pulsar in an 0.3 yr orbit this is already the case af- 
ter a few orbits. After a few years the Einstein delay can 
be measured with high precision, as will be shown in the 
simulations below. The dragging of inertial frames in the 
vicinity of the black hole also affects the Einstein delay. 
However, this occurs only at higher orders (Wex 1995), 
which in most cases can be neglected or easily accounted 
for in a combined mass and spin measurement. 

The Shapiro delay accounts for the extra light travel- 
ling time due to the curvature of space-time caused by the 
existence of surrounding masses (here mainly Sgr A*). 
The Shapiro delay contains two separately measurable 
PK parameters, the mass of the black hole Mbh and sin i. 
The signal is usually only sufficie ntly strong for edge-on 
systems (e.g. IKramer et al.l [20061 ), but in our case even 
for a face-on orbit (i = 0) the effect will be significant 
due to the large mass of Sg;r A*, if the orbit is eccen tric. 
Using the equation of iBlandford fc Teukolskvl (| 19761 ), 
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(7) 



as a first order estimation, one can see that for an ec- 
centricity of 0.5 the Shapiro delay for z = amounts 
to about 40 seconds. This already indicates that the 
Shapiro delay allows a precise mass determination, even 
for a pulsar with poor timing precision. Apart from con- 
taining Mbh directly, the Shapiro delay gives a second, 
though indirect, access to the Sgr A* mass via s'mi and 



^ The Roemer delay is defined as Ar = — Kq -r, where Kq is the 
unit vector of the hne-of-sight and r is the position vector of the 
pulsar with respect to the barycenter of the binary system. The 
Roemer delay describes the contribution of the pulsar motion to 
the time delay. 



the mass function. One finds 

where x is the projected semi-major axis of the pulsar 
orbit (in light seconds), which is an observable Keplerian 
parameter. It depends on the orbital eccentricity and 
inclination, which of the two is more constraining. 

In addition, there are significant contributions to the 
signal propagation caused by frame dragging. A first 
order analytic eq uation for this effect can be found in 
IWex fc KopeikinI (p^999). From this it is clear that the 
frame dragging can have a significant contribution to the 
propagation delay, but in most cases will have a distinct 
signature that can be fitted for, leading at the same time 
to a precise mass measurement and a lower limit on the 
spin parameter x- Contributions from higher order mul- 
tipole moments and light bending effects can easily be 
accounted for in an analytic way (see e.g. Kopeikin 1993). 

The inclination of the pulsar orbit with respect to the 
Hne-of-sight i (modulo a tt — z ambiguity, see Fig. S]) can 
be obtained either directly from the Shapiro delay, as ex- 
plained above, or via Eq. ([8]) by using the mass, Mbh, 
derived from any other PK parameter. Therefore, in 
Section [H and [5] where the determination of spin and 
quadrupole is presented, we can treat the inclination an- 
gle as a parameter that is known with sufficient precision. 
A brief discussion on the tt — i ambiguity can be found 
in Section l4Jl 

3.2. Simulations 

The simulations performed in this paper mainly con- 
tain two steps: creating TOAs and determining parame- 
ters as well as their measurement uncertainties. Firstly, 
the TOAs are created regularly regarding to solar sys- 
tem barycentric time and then combined with the three 
time delays (Roemer, Einstein and Shapiro, see the above 
Subsection) to account for the changes in the signal ar- 
rival time due to the orbital motion of the pulsar around 
Sgr A*. Next the simulated TOAs are passed to the 
TEMPO software package. Based on a timing model, 
TEMPO performs a least-square fit to yield a phase- 
connected solution of the TOAs, and determines the 
model parameters. The measurement uncertainties of 
these parameters are calculated via a covariance ma- 
trix. This is the standard procedure of pulsar timing 
observa tions, and explained in g reat detail in iTavlorl 
(I1994D: iLorimer fc Kramerl (|2QQ5 h: iHobbs eial] (j2QQ6[ ): 
lEdwards et al.l ( 20061 ). Most of the timing models used 
in this paper are part of the TEMPO standard imple- 
mentation available as a download from the sources given 
in these references. When ever we use an extension to 
these well tested models, to account for specific effects, 
which are not covered by the standard software, we will 
mention this explicitly in the corresponding Section. 

In this Subsection we present the simulations for the 
mass determination. For this we assumed five years of ob- 
servations with weekly TOAs which contain white Gaus- 
sian noise with a standard deviation of 100 /as. Fig. [3] 
shows the results of our simulations for a typical system 
configuration. If this is not the case then, as outlined 
above, u cannot a priori be used for a high precision mass 
measurement due to an unknown contribution from the 
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frame dragging, as we will show later. 

In practice, not only one single relativist ic effect will 
be used to determine the mass of Sgr A*, but a consistent 
model, accounting simultaneously for frame dragging ef- 
fects in the orbital motion and the signal propagation, 
will be used to determine the mass and spin at the best 
level. How the spin of Sgr A* affects the timing observa- 
tions and how it can be extracted from the timing data 
is the subject of the next Section. 
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Fig. 3. — Simulated fractional precision for the mass determi- 
nation of Sgr A* as a function of the orbital period Pb, obtained 
from three different relativistic effects: precession of the orbit (a;), 
Einstein delay (7e)5 and Shapiro delay (As). The mass determi- 
nations are based on simulated data, assuming weekly TOAs with 
an uncertainty of 100 /iS over a time span of five years. We used 
an orbital eccentricity e of 0.5 and an orbital inclination i, relevant 
for the Shapiro delay, of 60° . The simulations were done for a non- 
rotating black hole. Note that for various practical reasons (such as 
the uncertainty in the pulsar mass), a precision below 10~^ seems 
unrealistic. Also, as explained in the text, for a rotating black hole 
Cj cannot be used directly for a high-precision mass determination, 
due to the large contribution of frame dragging. 



4. FRAME DRAGGING, SPIN MEASUREMENT, AND GR'S 
COSMIC CENSORSHIP CONJECTURE 

Although there are clear indication that the Sgr A* 
rotates, its actual rate of rotation is still not well de- 
termined. Investigations of flares from accreting gas in 
the near-infr ared and in X-rays yield a ran^e of y ^ 
0.22 to 0.99 (i Genzel et al.' 2003; "Asche nbach et al.ll2004l : 
iBelanger et a l. 2006; Aschenbach 2010). The rather large 
range in the estimates of x is also a result of the uncer- 
tainty in the underlying model assumptions. A pulsar, 
however, would provide a clean probe of the gravitational 
field of Sgr A* and, in absence of any major external per- 
turbations, give a direct access to the dragging of iner- 
tial frames in the vicinity of Sgr A*. In iWex fc KopeikinI 
('1999) it has been shown that in relativistic pulsar-black 
hole binaries the (additional) precession of a pulsar orbit 
due to the frame dragging caused by the spin of the black 
hole (Lense-Thirring precession) is the most promising 
effect to determine the direction and magnitude of the 
black hole spin. This, in general, is also the case for a 



pulsar in orbit about S gr A*. The assumption made in 
IWex fc KopeikinI (|1999f ), that the spin of the black hole 
S is clearly smaller than the orbital angular momentum 
L, is no longer valid here. The fraction between the spin 
of the black hole and the angular orbital momentum is 
given by 

5'bh Mbh o /^^ 



where /3e = /3o/Vl — and MpsR is the mass of the 
pulsar. For a pulsar with an orbital period of less than 
one year is greater than 50000 x- Thus, the total 

angular momentum J of the system is completely domi- 
nated by the spin of Sgr A*, whose direction will therefore 
practically coincide with the direction of J, and can, for 
the considerations here, be viewed as a constant in time. 
In this case, the orbital motion to post-Newtonian accu- 
racy including first-order spin terms can be taken from 
Appendix B in We:^ (|1995D . This is sufficient to simulate 
all th e relevant effects (see iKonigsdorffer fc Gopakumarl 
l2005l for higher order corrections) for a system free of 
external perturbations. It also accounts for the fact that 
the precession is stronger near the pericenter. 

4.1. S'pin determination from the timing parameters 

Averaging over one orbit, one obtains the rates of the 
secular p recession of the pulsar orb it caused by frame 
dragging (jBarker fc 0'Connelll[T975[ n: 



$ = 
^ = 



— 3 l^LT COS^ 



Att 

VT = -^Plx^^X- (10) 



The definitions of the angles ^ and 6 are given in 
Fig. m The secular changes for the angles ^ and ^ are 
line ar in time. As d i scusse d in detail in IWexl (|1998f ) 
and IWex fc KopeikinI (|1999f ). this linear-in-time evolu- 
tion translates into a non-linear-in-time evolution of the 
observable angles that enter the timing model for a pul- 
sar, i.e. the longitude of pericenter uj and the inclination 
of the orbit with respect to the line-of-sight i (as part of 
the projected semi- major axis x). One findsl3 



and 



sin(cj — ^) 



cos{uj — ^) 



cx - ce Cj 
se Si 



(11) 



• (12) 



Since the angles z, 6>, and A are in the range to tt, 
their sines sx are non-negative and can be expressed as 
sx = \/l — Cx. As shown by Wex (1998), if the change 
in <l> is small (less than a few degrees) over the time span 
of the timing observations, the most straightforward way 
to analyze the timing data is to fit for the coefficients of 
the Taylor expansion of uo{t) and x{t) 

U=UQ^UQ(t-tQ)^\uQ(t-tQf ^ ... , (13) 

x = xo+xo(^-^o) + |^o(^-^o)^ + ... , (14) 

^ To estimate the measurability of the Lense-Thirring effect, it 
is sufficient to use the averaged precession rate. In practice, the 
preccession of the orbital plane is more complicated as can be seen 
from the analytic solution given in Appendix B of i Wexi (^1995i ). 

^ We define cx = cosX and sx = sinX. 
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and to use the parameters cjq, xq and their time deriva- 
tives as intermediate parameters to determine the angles 
6>, A, and the spin parameter x- For the config- 

urations considered in this paper it is sufficient to keep 
only terms up to second order in t — to. Nevertheless, 
we have extended TEMPO to account for cubic terms 
in oder to test their significance in all our simulations. 
We would like to note in passing, that the coefficients of 
the cubic terms can be calculated from the other coeffi- 
cients based on basic geometric relations, and therefore 
they would not add further information concerning the 
orientation of the system and the spin magnitude. 

From the derivatives of Eqs. (pT]) and (p!2]) one finds the 
relation between the time derivatives, the orientation of 
the orbit at a given epoch, and the spin of Sgr A*. In 
practice, the linear trend becomes visible in the timing 
data soon after the beginning time of observation, al- 
lowing the measurement of xq and the extraction of the 
Lense-Thirring contribution from c^o- One finds (for con- 
venience we drop the index 0): 



-XS- QSal^LT 



[(1 - 3«f 



)ce - CiCx 



LT 



(15) 
(16) 



where 53 = sqSxs^. Since at this point x, and Si 
are known quantities, the measurement of x determines 
the quantity \ss\x which must not exceed unity, since 
according to the CCC X ^ 1 ksl ^ 1 by definition. 
Hence, as soon as x becomes measurable one has a first 
test for the CCC. 

To fully determine the magnitude and orientation of 
the spin, the measurement of higher order derivatives is 
necessary. The second time derivatives read 



2 

LT 5 



x = -xs. ^ [si + s'^Ci{c0Cx - Ci)] Q 
uj = s~^ [2ciCe - (2 - 5-)ca] ssI^lt ^ 



(17) 
(18) 



which give us now, in total, six equations (pT| [12) [T5l 
[Ip [la [m for five unknowns (l9. A, ^0, ^0, x)- For a 
discussion of the solution of these equations, we introduce 
the variables xe = cex, Xx = caX, and (s = ssX- The 
parameters xe and xa represent the projection of the 
(normalized) spin onto the orbital angular momentum 
and the line-of-sight direction, respectively. They can be 
determined from the timing parameters via the equations 
of the first time derivatives 



XS7 



{uj - ujm )s: 



xQ 
2 



(19) 



m = {l-3s^^)xe-c^xx. (20) 



and those of the second time derivatives 
{xxc^ x'^s'^)sj _ 



{xttciy 



^ ^2 = q {Xe + Xa) - Ci(l + )x^Xa , 



xCt 



W2=1c\xe-Ci(\^c\)xx 



(21) 
(22) 



where ^3 has been eliminated using Eq. (jT9|) . The quan- 
tities X\ , Wi , X2 , and W2 are defined such that they do 
not change when the sign of Ci is fiipped. The above equa- 
tions can be easily solved analytically. By the time the 



second derivatives are observable, s^, i, cL^lt, and so X\ 
and >Vi, will be known with high precision. For a given 



sign of Ci = ±Y^1 — 5^, Eq. (j22j) will lead to a unique so- 
lution for xe^ Xa, and (s- For some orientations Eq. (|2T]) 
turns out to be more constraining. However, this gives us 
in general two solutions for (x6'7Xa)- But then Eq. (|22]) . 
although less constraining, can be used to rule out one of 
the two solutions. The best way to represent the solution 
is to plot the constraints from Eqs. (|2n]) . (|2T]) and (|22]) in 
the X9~Xx plane. Possible solutions are represented by 
the region where all three curves meet within the uncer- 
tainty given by the measurement errors of i, cL^lt, ^, and 
Co. This will become clear in the Section below, where 
we present the simulations. 

With xoi Xx and (s known, we can calculate the spin 
parameter of the black hole via 



X = yJCi^Xe^xl- '^CiXeXx • 



(23) 



Once X is determined, we can calculate all the angles. Fi- 
nally, the i ^ TT — i ambiguity leaves us with two different 
solutions in the orientation by (<l>o. A) ^ (tt + <l>o, tt — A), 
for which x has the same value. 




pericenter 



Fig. 4. — Definition of angles in Sgr A* spin reference frame. The 
orientation of the orbit with respect to the observer is given by the 
orbital inclination i and the longitude of pericenter uu as measured 
from the ascending node in the plane of the sky. The pulsar orbit 
with respect to the equatorial plane of the rotating black hole is 
determined by the inclination 0, the equatorial longitude of the 
ascending node <l>, and the equatorial longitude of pericenter ^. 
The angle between the line-of-sight and the Sgr A* spin is denoted 
by A. 



4.2. Simulations 

The technique described in the previous Subsection has 
been tested by a set of standard simulations for various 
orbital conffgurations. For a given system, following the 
procedures described in Section 13.21 we simulate weekly 
100 /is TO As over a time span of ffve years. Here in the 
calculations of the time delays, in addition to the rela- 
tivistic pericenter advance we also consider the inffuence 
of Sgr A* spin by inputting the secular changes of $ and 
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^ described in Eq. ([TO]). Then, in order to determine the 
PK parameters, we fit the TOAs with the MSS timing 
model of TEMPO, which we have extended to model the 
secular changes in pericenter and projected semi-major 
axis up to the third order in the time derivatives. The 
third order coefficient turn out not to be significant in 
the simulations presented here. Figs. [5] and [6] show the 
Xe~Xx plane for two different orientations of the black 
hole and the pulsar orbit. According to GR the solu- 
tion has to lie within the boundaries of the figures, since 
— 1 < X0,Xa < 1 for a Kerr black hole. Moreover, the 
solution (x6>,Xa) has to lie within an ellipse defined by 
setting X = 1 in Eq. (|23|) . in order to represent a Kerr 
black hole with an event horizon. Once x is measured, 
one can determine ("s from Eq. (p!9|) and use this quantity 
to plot the ellipse defined by Eq. ([23]) in the X6>-Xa plane. 

Fig. [71 shows a simulation for a Kerr solution with a 
spin that exceeds the spin of an extreme Kerr black hole. 
Within GR, this would represent a naked singularity. For 
such an object the CCC is violated and the predictabil- 
ity of the (classical) theory breaks down. Also, all three 
lines have to agree in a common region, otherwise either 
GR is not the correct theory, or there are external per- 
turbations present, a situation which we discuss in more 
details in the next Section. 




-1 -0.5 0.5 1 



Xe = X cos(e) 

Fig. 6. Like Fig. [S] but $o = 105°, 6> = 30°, and A = 75°. 
The corresponding spin parameter, as calculated from Eq. (|23|) , is 
X = 1.0001 ± 0.0003 (95% C.L.). 
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Fig. 5. — Determination of the SgrA* orientation in the xe~ 
Xx plane. For this simulation we have used an orbital period of 
0.3 yr, an orbital eccentricity of 0.5, x = 1, = 45°, <l>o = 45°, 
= 60° , and A = 60°, which are in agreement with the constraints 
by Zamaninasa b et al.l (|201in . A change in the sign of q mirrors 
the figure along the XA = line, meaning that the solution for 
is invariant, but A changes to tt — A. The corresponding spin 
parameter, as calculated from Eq. 03, is X = 0.9997±0.0010 (95% 
C.L.). In all the X6'~Xa plots (Fig.[5]-[8l) we plot the 68% confidence 
intervals. However, in most cases the separation between the two 
lines is below the resolution of the plot. The dotted ellipse is the 
boundary of the area for Kerr black holes (see text for details). 



4.3. Identification of external perturbations 



Fig. 7. — Parameters as in Fig. \E\ but x — 1-2 (naked Kerr 
singularity). The dotted ellipse is the (outer) border of the region 
where, for the measured orbital inclination and x, the Kerr black 
holes are located, i.e. where x ^ 1- 

As discussed in great detail by iMerritt et al.l (|2010[ ), 
the orbit of a star or pulsar around SgrA* may be sub- 
ject to perturbations from other stars in the vicinity of 
the black hole. Depending on the number density of the 
stars, this could significantly affect the precession of the 
pulsar orbit. Nevertheless, since we have three lines in 
the X9~Xx plane that need to intersect, our analysis will 
unveil the presence of any external perturbations. In 
Fig. [8] we present a X9~Xx diagram based on timing data 
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that contain (besides the gravitational field of Sgr A*) 
an external perturbation causing an additional preces- 
sion of the pericenter. For orbits with Pb ^ 0.3 yr, even 
a small (compared to the Lense-Thirring precession) ex- 
ternal contribution to the precession of the pericenter 
leads to a situation where the cj, C(j, and x lines fail to 
intersect in one point within the measurement precision. 
The same is true, if there is an external contribution to 
a change in the inclination of the orbital plane. Hence, 
if all three lines intersect, we not only have a precise de- 
termination of the spin of the black hole, but also a test 
that this measurement is not contaminated by external 
perturbations. 

At this point we would like to add a more detailed com- 
ment on the discriminating power of the pulsar test con- 
cerning external perturbations. In practice, the three- 
line test outlined above is not simply based on the sec- 
ular precession rates. We emphasize that a consistent 
fit of the timing data, with a model that includes the 
Lense-Thirring precession, needs to incorporate the full 
dyna mics o f the orbital precession as given by Appendix 
B in IWexl ([19951). The phase dependence of the Lense- 
Thirring precession rate is a direct result of the Coriolis 
type force caused by the dipolar gravitomagnetic field of 
the central rotating black hole. Hence, we can identify an 
external perturbation based on this quasi-periodic effect, 
even in a fine tuned situation where the external mass 
distribution manages to mimic a secular Lense-Thirring 
precession. In fact, we have conducted simulations and 
found that the phase dependent precession rate leads to 
an effect that is typically four orders of magnitude larger 
than the timing precision assumed in our simulations. 
This is in line with the findings of iDamour fc Deruelld 
(|l986l ). who pointed out the strength of quasi-periodic 
effects in tests of gravity. 
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Fig. 8. — Parameters as in Fig.O but the precession of pericenter 
has an additional contribution from an external perturbation that 
amounts to 10% of the Lense-Thirring contribution. For a better 
resolution only the first quadrant of Fig. [5] is plotted here. 



5. QUADRUPOLE MEASUREMENT AND NO-HAIR 
THEOREM TEST 

The quadrupole moment of the black hole leads to an 
additional secular precession of the pulsar orbit. This 
precession, however, is even for compact orbits (P5 ~ 
0.1 yr) much smaller than the Lense-Thirring preces- 
sion. Further, it can be shown that the secular terms 
of the precession cannot be separated from the Lense- 
Thir ring effect. Fo r this reason, it has been argued 
by iWex fc Kop eikin (1999) that while the spin magni- 
tude and the orientation of the black hole are mainly 
determined by the overall precession of the orbit, the 
quadrupole of the black hole is mostly determined via its 
periodic infiuence on the motion of the pulsar from one 
pericenter to the next. As will be shown in this Section, 
these periodic features of the quadrupole can be used to 
fit for the quadrupole moment of Sgr A*. 

5.1. Extracting the quadrupole from the timing data 

The deviations in the motion of the pulsar caused by 
the quadrupole moment lead to a variation in the Roemer 
delay, which we describe by a change in the coordinate 
position of the pulsar according to 

= (r + Jr(^))(n + Jn^^)) . (24) 
The vector is calculated from the changes in the angles 

= <!> + (5<l>(^) , = <!> + (5^^^) , e' = e^ 60^"^^ , 

(25) 

according to Sh = — n. To first order in e = 
— 3Q/a^(l — e^)^, the detailed equations for the ^- 
quantities can be taken from Garfinkel (1959), with slight 
modifications that account for the dominating precession 
of pericenter caused by the mass monopole: the term 
(5?/^ — 1) in the auxiliary constants m and 7 has to be re- 
placed by 2cjPb/7re, where uj is the total advance of peri- 
center. Based on this we have developed a timing model 
that includes the contribution of the quadrupole moment 
of Sgr A* to first order in e. Fig. [9] illustrates the unique 
periodic timing residuals caused by the quadrupole mo- 
ment of Sgr A*. 




Fig. 9. — Residuals caused by the quadrupole moment of Sgr A* 
plotted for two orbital phases. We have used the same orbital and 
black hole parameters as in Fig. \E\ 
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This periodic signal will not only allow the determi- 
nation of the quadrupole moment of SgrA* with high 
precision, but also provide a clear identification of the 
quadrupolar nature of the gravitational field. Moreover, 
due to the large advance of pericenter the quadrupolar 
signal will change in a characteristic way from one orbit 
to the next. This clearly helps to identify any external 
"contamination" of the orbital motion of the pulsar, and, 
as in the spin determination, provides high confidence in 
the reliability of a no-hair theorem test with a pulsar 
around SgrA*. 

5.2. Simulations 

We have tested the procedure outlined above in a num- 
ber of mock data simulations, for various orbital config- 
urations^ Again following the procedures described in 

Section 13.21 we assume weekly TOAs with a precision of 
100 /is for a time span of five years. This time we ex- 
tended our simulations and the timing model used in 
Section 14.21 to account for the periodic effects due to 
quadrupole moment of SgrA* described in Eq. Our 
results are summarized in Fig. [TOl Note that the preci- 
sion of the spin determination is at least one order of 
magnitude better than the determination of q. Hence 
the uncertainty in the g'-measurement is the limiting fac- 
tor for the no-hair theorem test. As a conclusion of our 
simulations, if the external perturbations are negligible, 
for orbits with Pb ^ 0.5 yr the no-hair theorem can be 
tested with high precision. If we adopt the precessional 
rates from stellar perturbation calculated in Fig. [21 we 
conclude that the test can be achieved with high pre- 
cision for orbits with Pb ^ 0.1 yr. This range can be 
extended if in the presence of perturbations the char- 
acteristic quadrupolar features remain separable. This, 
however depends on the details of the external mass dis- 
tribution, which we will not investigate further in this 
paper. 

6. DISCUSSION 

In this paper we have developed a method to deter- 
mine the mass, the spin, and the quadrupole moment 
of Sgr A* using a pulsar in a compact orbit around this 
super-massive black hole. Our investigation is based on a 
consistent timing model, that includes all the relativistic 
and precessional effects that can be used to extract these 
parameters of SgrA*. Based on simulated timing data 
for a pulsar in orbit around SgrA*, we have shown in 
a consistent covariance analysis, that even with a mod- 
erate timing precision (~ 100 /is), one can expect to be 
able to determine the mass, the spin, and the quadrupole 
moment of SgrA* with high precision, provided the or- 
bital period is well below one year. As a result of our 
simulations, for a compact orbit (orbital period of a few 
months) one can expect to measure the spin with a pre- 
cision of 10~^, or even better. We have shown how the 
method would allow the identification of an object whose 
frame-dragging exceeds that of an extreme Kerr black 
hole, and therefore would provide a test of the CCC. 
Furthermore, for such orbits the determination of the 
quadrupole moment of SgrA* seems feasible at a few 
percent precision level or even better, depending on the 
size and orientation of the pulsar orbit and the spin of 
SgrA*. In combination with the precise spin measure- 
ment from the Lense-Thirring effect, this yields a high 




Pb [yr] 

Fig. 10. — Measurement precision for the quadrupole moment of 
Sgr A* as a function of orbital period for three different eccentric- 
ities, in absence of any external perturbations. We have used the 
same orbital and black hole parameters as in Fig. [5] For the tim- 
ing, we assumed the same time span and characteristics of TOAs 
as in Fig. O This time however the TOAs were equally distributed 
with respect to the true anomaly, in order to account for the fact 
that timing needs to be done more frequently around pericenter to 
optimize the measurement of the quadrupolar signal in the TOAs. 



precision test of the no-hair theorem of stationary black 
holes. 

Moreover, we have shown that, in general, our analysis 
will be able to unveil the presence of external perturba- 
tions caused by the presence of distributed mass, there- 
fore providing high confidence in a spin and quadrupole 
determination based on pulsar timing. Nevertheless, fur- 
ther studies are required to see whether the spin and 
quadrupole moment can still be extracted if the tim- 
ing data is "contaminated" by external perturbations. If 
perturbations arise from a smooth concentration of dark 
matter particles in the vicinity of SgrA*, we may be able 
to learn something about the properties of dark matter 
that clusters around SgrA*, assuming GR is correct. 

Finally, we need to emphasize that the tests presented 
are not affected by an uncertainty in the distance to the 
GC. On the contrary, a mass determination via pulsar 
timing would give a greatly improved value for Rq if 
combined with the astrometric measurements in the near 
infrared. 

Once a pulsar is detected in a compact orbit around 
Sgr A*, continuous timing will allow more and more mea- 
surements and tests as the timing baseline grows with 
time. In the following we summarize the most important 
steps in this experiment: 

• After timing one orbit, all Keplerian parameters 
will be well known, and also the pericenter ad- 
vance will be measured with good precision. This 
will already provide a good estimate of the mass of 
SgrA*. 

• Timing a few more orbits would then allow the de- 
termination of additional post-Keplerian parame- 
ters, like the Shapiro parameters (rsh, sinz) and 
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the ampHtude of the Einstein delay (7e)- These 
parameters allow a robust determination of Sgr A* 
mass and the inclination of the pulsar orbit with 
respect to the Hne-of-sight. 

• At this stage, a measurement of a change in the 
projected semi-major axis (x) caused by the Lense- 
Thirring will allow an early test of the CCC, and 
mapping of the region in the X9 ~ Xx plane where 
the solutions for Kerr black holes are. 

• Around the same time the mass measurement 
should reach a precision that allows the extraction 
of the Lense-Thirring contribution to the preces- 
sion of the pericenter (cOlt), giving a line-like re- 
gion in the xe — Xx plane. 

• After a few years of timing the second time deriva- 
tives of u and x should be known with high preci- 
sion, allowing a precise determination of the Sgr A* 
spin (magnitude and direction). At this stage we 
also have a test for the "cleanness" of the sys- 
tem, and whether the spin is below the Kerr bound 

(x = 1). 

• At the same time the obtained parameters for the 
pulsar orbit and the Sgr A* spin can be used to 
model the periodic features in the timing residu- 
als, which are caused by the quadrupole moment 
of Sgr A*. This allows to a determination of the 
quadrupole moment and a test of the no-hair the- 
orem. 

A potential problem for the timing of a pulsar in 
a compact orbit around Sgr A* is posed by the rela- 
tivistic spin precession, as pointed out by iMerritt et al.l 
()2Q1G) . This change in the pulsar orientation with re- 
spect to a distant observer not only causes a variation of 
the pulse profile, which makes precise timing more diffi- 
cult, but also ca n turn the pulsa r emission away from 
our line-of-sight (jWeisberg et al.l 119891 : iKramerl 11998). 
To leading order the spin-precession is given by the de 
Sitter precession rate, which for Mbh ^ ^psr reads 
([Barker fc O^Connelllll975[ ) 



dS 



/3l ^ (0.137yr) 



lyr 



-5/3 



(26) 



Consequently, for orbital periods below one year rela- 
tivistic spin precession is expected to play an important 
role in the timing observations. We note in passing, that 
the Pugh-Schiff precession rate caused by frame drag- 
ging (l^FD - 2^/?gv/PK. lPughlll959l : lSchiffl ll960) is only 
relevant in the case of very compact, highly eccentric or- 
bits (^^FD ^ I7yr for Pb = 0.1 yr, e = 0.8 and x = 1), 
and could provide an independent test of the rotation of 
Sgr A*. 

The no-hair theorem test can also be affected by the ac- 
cretion disc around Sgr A*. To get an idea of the strength 
of this effect, one can estimate the influence by calculat- 
ing the fraction of the quadrupolar potential of the disc 
to that of the black hole. This ratio is given by 



MdiscrLc/(^BHr| 



(27) 



where Mdisc and raise are the mass and outer radius of 
the disc, and Vg = GMbu/c^ is the black hole gravita- 
tional radius. Following the advection -dominated accre- 
tion flow (ADAF) model of .Yuan et al.l (|2QQ9) and adopt- 
ing, as an upper limit, the disc scale of ^ 1 arcsec d eter- 
mined from X-ray observation (Baganoff et al.l l2QQlf ). we 
obtain ^ ^ 0.4%, which indicates that the quadrupole 
moment measurement of Sgr A* would not be biased by 
the contribution of the disc above the 1% precision level. 
Furthermore, in a very resent publication, based oncur- 
rent X-ray and millimeter observations, iPsaltisI ()2Qlll ) 
concludes that for compact orbits, like the ones discussed 
in this paper, hydrodynamic drag forces from plasma in 
the vicinity of Sgr A* are expected to be negligible. 
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